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Abstract. The generation of large scale flows by the anisotropic kinetic alpha (AKA) effect is investigated in 
simulations with a suitable time-dependent space- and time-periodic anisotropic forcing lacking parity invariance. 
The forcing pattern moves relative to the fluid, which leads to a breaking of the Galilean invariance as required for 
the AKA effect to exist. The AKA effect is found to produce a clear large scale flow pattern when the Reynolds 
number, 1Z, is small as only a few modes are excited in linear theory. In this case the non- vanishing components of 
the AKA tensor are dynamically independent of the Reynolds number. For larger values of 7Z, many more modes 
are excited and the components of the AKA tensor are found to decrease rapidly with increasing value of 1Z. 
However, once there is a magnetic field (imposed and of sufficient strength, or dynamo-generated and saturated) 
the field begins to suppress the AKA effect, regardless of the value of 1Z. It is argued that the AKA effect is 
unlikely to be astrophysically significant unless the magnetic field is weak and 1Z is small. 
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1. Introduction 

Vorticity and magnetic fields display some important sim- 
ilarities. Both satisfy formally similar equations: the vor- 
ticity and induction equations, respectively. This analogy 
has been used extensively by Batchelor (1950) in early 
work on hydromagnetic turbulence and dynamos. Indeed, 
both vorticity and magnetic field vectors are found to 
be aligned with each other in simulations of convection 
with dynamo-generated magnetic fields (Brandenburg et 
al. 1996). However, the intensity of the two vector fields 
does not generally show any correlation, and even slightly 
different initial conditions for vorticity and magnetic field 
lead to their mutual departure after some time (Moffatt 
1978). 

A formal analogy between vorticity and magnetic fields 
has also been proposed in the context of mean-field the- 
ory. Moiseev et al. (1983) proposed the possibility of an 
alpha effect in the equation for the mean vorticity in a 
compressible fluid. However, if a similar effect is to exist 
in the incompressible case there must be some anisotropic 
forcing, because otherwise the non-dissipative terms in the 
equation for the mean velocity would vanish (Krause & 
Riidiger 1974, see also Moffatt & Tsinober 1992). 

Frisch et al. (1987) were the first to study the effects of 
a non-Galilean invariant forcing that produced a destabi- 
lization of the velocity field at large scales. In analogy with 
the alpha effect in dynamo theory (Moffatt 1978) they 
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called this the anisotropic kinetic alpha effect, or AKA ef- 
fect. The nonlinear behaviour of this effect was studied 
by Sulem et al. (1989) and Galanti & Sulem (1991), who 
showed that the energy transfer to larger scales happens 
successively, in the form of an inverse cascade. 

In recent years a number of geophysical and astrophys- 
ical applications of the AKA effect have been discussed in 
the literature. For example, Zimin et al. (1989), Rutkevich 
(1993), and Levina et al. (2000) applied the effect to geo- 
physical convection. Krishan (1991) discussed the possi- 
bility of an inverse cascade in the context of solar granu- 
lation and Krishan (1993) invoked the anisotropic kinetic 
alpha effect in order to explain the clustering of galaxies. 
Kitchatinov et al. (1994) and v. Rekowski et al. (1999) 
discussed the possibility of creating large scale vortices in 
accretion discs. Recently, v. Rekowski & Riidiger (1998) 
suggested that the AKA effect could help to solve the 
'Taylor number puzzle' in models of stellar differential ro- 
tation. They found that this effect could produce angular 
velocity contours that are no longer constant on cylinders, 
and hence closer to the helioseismological observations. 

In magnetohydrodynamics the generation of large scale 
magnetic fields has been well established numerically in 
simulations of the geodynamo (Glatzmaier & Roberts 
1995) and accretion discs (Brandenburg et al. 1995), for 
example. However, the excitation condition for the gener- 
ation of large scale vortices should be similar to the con- 
dition for dynamo action (Kitchatinov et al. 1994). It is 
therefore remarkable that no evidence for the spontaneous 
generation of vortices has been seen in the simulations 
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of Brandenburg et al. (1995), even though the boundary 
conditions for vorticity and magnetic field were identical. 
The same is true of the recent simulations of Brandenburg 
(2001) using fully helical isotropic flows in periodic do- 
mains. No trends for large scale flows are seen, even though 
the magnetic field displays a very pronounced large scale 
pattern. 

The evidence for large scale flows is sparse. However, 
one example might be the large scale flows seen in 
highly supercritical Rayleigh-Benard convection (Howard 
& Krishnamurti 1986). Here it could be the presence of 
boundaries which breaks Galilean invariance. By contrast, 
the simulations of Brandenburg (2001) were Galilean 
invariant and isotropic, which explains the absence of 
an AKA effect. Astrophysical examples where turbu- 
lence is driven by non-Galilean invariant forcing include 
supernova-driven turbulence in galaxies (Korpi et al. 
1999) and the turbulent wakes driven by galaxies moving 
through the galaxy cluster (Ruzmaikin et al. 1989). 

It is important to realize that the AKA effect has been 
verified numerically only in the case of rather low Reynolds 
numbers, 1Z < 2. Given the possible astrophysical rele- 
vance of the AKA effect, it is important to assess the de- 
pendence of the resulting large scale flow on the Reynolds 
number, which is extremely high in all astrophysical set- 
tings. Another important property of astrophysical flows 
is that the gas is ionized and electrically conducting, so it 
may be unstable to dynamo action. Typically, the result- 
ing large scale magnetic fields attain near-equipartition 
strength and will therefore be dynamically important. The 
combined action of AKA and a effects has already been 
considered by Galanti et al. (1990, 1991), but again only at 
relatively small kinetic and magnetic Reynolds numbers. 

The purpose of the present paper is to extend the stud- 
ies of Frisch et al. (1987) and Sulem et al. (1989) to the 
case of larger Reynolds numbers, allowing also for a mag- 
netic field to grow from a weak initial seed magnetic field. 
In all cases we adopt the same forcing as Frisch et al. 
(1987). However, in contrast to their original paper, where 
the flow was assumed to be incompressible, we assume 
here weak compressibility. 

2. Description of the model 

We solve the isothermal compressible hydromagnetic 
equations for the logarithmic density In p, the velocity u, 
and the magnetic vector potential A, 



D hip 
Dt 

Dm 

TJT _ 

dA 
~dt 



V • u. 



(1) 



- c 2vinp+^^ + ^(V 2 t t +iVV-t t ) + /, (2) 



u x B — rjfioJ, 



(3) 



box, D/Di = d/dt + u ■ V the advective derivative, B = 
V x A the magnetic flux density, J — V x B/po the 
current density, c s the sound speed, and 



/ = V2./o <f> 2 (4) 




is the forcing term of Frisch et al. (1987) with 

4>i = cos(fcfj/ + u>t), <f>2 = cos(k{X — tot), (5) 

and u> — vkf. This forcing corresponds to a pattern mov- 
ing with the velocity (1, —1, 0)vk{ diagonally in the (x, y)- 
plane. k{ is the wavenumber of the small scale forcing and 
fo (= const) gives the strength of the forcing. The (un- 
curled) induction equation (||) implies a specific gauge for 
A such that the electrostatic potential vanishes. Instead 
of the dynamical viscosity p (— const - not to be confused 
with the magnetic permeability po) we will in the follow- 
ing refer to the mean kinematic viscosity v = p/po, where 
Po is the volume averaged density (po = const owing to 
mass conservation). 77 (= const) is the magnetic diffusivity. 

In order to nondimensionalize the equations, velocity 
is measured in units of vkf, length in units of kf , time 
in units of {vkf)^ 1 (= lj^ 1 ) and magnetic field in units 
of ^/poPo^fcf • The nondimensional form of Eq. ([!]) is then 
unchanged, and the nondimensional momentum and in- 
duction equations become 



K 2 

— = -w vlnp 



Dm 

dT 



J x B V 2 « + ±VV-u ^ „ f . 
+ +ft/,(6) 



dA 1 T 

— — = u / B —-J- 

3t V m 



(7) 



in a three-dimensional periodic cartesian domain of size 
L = 2n/ki, where k\ is the smallest wavenumber in the 



where all variables are nondimensional, and the nondi- 
mensional forcing function is / = \/2 (0i, 4>2,4>i + 1^2) ■ I 11 
the following, however, we express all relevant variables in 
explicitly nondimensional form, e.g. we write cot instead 
of just t. 

The problem is completely defined by four nondimen- 
sional parameters: the Reynolds number 1Z = Uo/{vkf), 
where Ua = Jq/lo is a reference velocity, the magnetic 
Reynolds number lZ m = Uo/(rjk{), the Mach number 
A4 = Uq/c s , and the size of the computational domain 
L, which we express in terms of the nondimensional quan- 
tity fcf/fci = k{L/(2ir). The magnetic Prandtl number is 
V m = !Z m /lZ. In some additional runs we also applied an 
external field Bo — Boz, which leads to a fifth parameter 
Bq = B /(^/fi p Uo)- We note that our definition of 1Z 
agrees with that of Frisch et al. (1987). 

As a representative measure of the resulting veloc- 
ity field we monitor the normalized rms velocity, hi = 
(it 2 ) 1 / 2 /Uo- Here the angular brackets denote averaging 
over the volume of the domain. The magnetic field is mon- 
itored analogously through B — (B 2 ) 1 / 2 / (^/poPo^o) +&o, 
where Bq quantifies the imposed magnetic field. However, 
in most of the cases where we have a magnetic field we rely 
on the dynamo-generated field and put Bq = 0. We also 
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calculate the normalized rms velocity of the large scale 
flow, Uhs = (u 2 ) 1 ' 2 /Uq. With the forcing function chosen 
the z-direction is preferred and horizontal (x, y) averages 
(denoted by overbars) are appropriate for extracting the 
large scale how. We then give the ratio Uls/U. Wc note, 
however, that a definition of Wls in terms of Fourier fil- 
tering to k = k\ is also sensible. In some cases we give the 
relative kinetic helicity of the horizontally averaged (large 
scale) flow, 



o_/LS 

n K 



u)/{{tf){u 2 )f' 2 , 



(8) 



where angular brackets indicate volume averages (as op- 
posed to the overbars which denote only horizontal aver- 
ages). We also calculate the nondimensional growth rate 
of the magnetic field, \/u>, the ratio of rms magnetic field 
to rms velocity, B /U, as well as nondimensional measures 
of kinetic and magnetic hclicities, 



Hk = (u-u)/((u> 2 )(u 2 )) 1 / 2 1 
Hm = {A-B)/({A 2 c ){B 2 )Y/ 2 , 



(9) 
(10) 



respectively. Here, A c = A — A$ is the magnetic vector 
potential in Coulomb gauge, where A$ = (A) + V</> and 
V 2 = V • A. The advantage of using A c is that it has the 
property of minimizing (A 2 ). The quantity (A-B) is gauge 
invariant in periodic domains (and therefore (A c ■ B) = 
(A ■ B), for example). If there is an additional imposed 
uniform magnetic field, (A ■ B) would no longer be gauge 
invariant, because (A) ^ in general. Since a uniform field 
in isolation has zero magnetic helicity we calculate (A ■ B) 
only with respect to the deviations from the imposed field. 
Finally, we also monitor the quantity 



H M = k 1 (A-B)/(B 2 



(11) 



which is a nondimensional length scale. The reason for 
using here k\ instead of ki is that the magnetic helicity 
tends to develop on the largest possible scale, independent 
of the forcing wavenumber. 
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Fig. 1. Upper panel: Power spectra of kinetic energy, 
E(k,t), for Run 9c at three different times. Note the 
gradual build-up of power in the lowest wavenumber, 
k\ = fcf/14. For reference, the dissipative cutoff wavenum- 
ber, defined as kd = {{oj 2 ) /v 2 ) 1 ^ , is indicated by arrows 
in the upper- right corner separately for the three times. 
Lower panel: E(k) (at t — 660/w) is compared with the 
helicity spectrum, F(k), normalized by 2k. The diamonds 
mark the few points where F(k) > 0; everywhere else 
F(k) < 0. Note that only at the smallest wavenumber the 
helicity is close to its largest possible value < 2kE by 
the realizability condition). 7Z — 8. 



3. Results 

We first consider the case where TZ m is so small there is 
no dynamo action. We have calculated solutions for differ- 
ent values of TZ (between w 0.82 and 12) and kf/k\ 
(between 6 and 14). The case TZ = ^/2/3 and fcf/fci = 6 
was considered by Frisch et al. (1987) for the strictly in- 
compressible case and neglecting any magnetic fields com- 
pletely. Without magnetic fields of sufficient strength we 
find in all cases clear signs of the AKA effect, which man- 
ifests itself through enhanced power in the lowest Fourier 
modes of the kinetic energy spectrum; see Fig. |l|. It turns 
out that at least for small values of TZ the spectral power 
in the lowest Fourier mode exceeds that of the forcing 
wavenumber. As the value of TZ is increased the solution 
becomes more irregular in time (Fig. ||), as also found by 
Sulem et al. (1989). It is also noteworthy that, at least 
for moderately large values of TZ, the spectral peak at the 



smallest wavenumber is relatively broad compared with 
the magnetic inverse cascade where, again using a peri- 
odic box, the peak is very sharp; cf. Figs. 17 and 19 of 
Brandenburg (2001). 

A summary of all the runs that we have performed is 
given in Table [l|. Here we list some characteristic quan- 
tities as a function of TZ and other input parameters. 
Perhaps the most important diagnostic quantity is the ra- 
tio Uls/U characterizing the relative importance of the 
resulting large scale flow. As a rule, when Uls/U > 0.7, 
the spectral energy at the largest scale exceeds that at 
the forcing scale. For U^sjU in the range 0.4 to 0.5 the 
two are comparable on average, but the large scale flow is 
unsteady in time. For U\sjU less than about 0.3 there is 
usually only the one peak at the forcing scale. This is the 
case especially when there is a magnetic field, i.e. B/U is 
finite. 
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Table 1. Summary of the main properties of various runs. A hyphen in the entry for B indicates that the magnetic 
field decays. 



Run 


mesh 


kf/ki 


TZ 


n m 


M 


Bo 


u 


Wls/W 


X/lo 


B/U 


Tin 




T~Cm 


la 


32 3 


6 


0.82 


0.82 


0.15 





2.40 


0.86 


-0.04 


— 


-0.19 


— 


— 


lb 


32 3 


6 


0.82 


41 


0.15 





1.36 


0.09 


0.05 


0.24 


-0.06 


0.00 


0.00 


lc 


32 3 


6 


0.82 


0.82 


0.15 


0.07 


2.10 


0.78 


0.02 


0.24 


-0.16 


— 


— 


Id 


32 3 


6 


0.82 


0.82 


0.15 


0.14 


1.67 


0.60 


0.01 


0.39 


-0.09 


— 


— 


2a 


30 3 


6 


1 


0.33 


0.10 





1.46 


0.73 


- 


- 


-0.14 


- 


- 


3a 


30 3 


6 


1.2 


0.4 


0.12 





1.60 


0.81 


- 


- 


-0.19 


- 


- 


4a 


30 3 


6 


1.5 


0.5 


0.15 





1.61 


0.85 


- 


- 


-0.21 


- 


- 


4b 


30 3 


9 


1.5 


1.5 


0.15 





2.91 


0.87 


-0.02 


- 


-0.26 


- 


- 


4c 


30 3 


9 


1.5 


15 


0.15 





1.55 


0.08 


0.15 


0.84 


-0.04 


0.35 


0.19 


5a 


30 3 


6 


2 


0.67 


0.20 





1.52 


0.88 


- 


- 


-0.22 


- 


- 


5b 


60 3 


14 


2 


2 


0.18 





1.77 


0.44 


-0.01 


— 


-0.08 


— 


— 


5c 


60 3 


14 


2 


20 


0.18 





1.29 


0.04 


0.20 


0.83 


-0.05 


0.29 


0.11 


6a 


30 3 


6 


2.3 


0.77 


0.23 





1.48 


0.89 






-0.23 






7a 


30 3 


6 


3 


1 


0.30 





1.43 


0.84 


- 


- 


-0.23 


- 


- 


7b 


60 3 


6 


3 


1 


0.06 





1.78 


0.77 


0.00 


- 


-0.20 


- 


— 


7c 


60 3 


9 


3 


1 


0.06 





1.43 


0.41 


-0.06 


- 


-0.10 


- 


- 


7d 


60 3 


14 


3 


1 


0.06 





1.56 


0.32 


-0.03 


- 


-0.09 


- 


- 


7e 


60 3 


9 


3 


- 


0.17 





1.60 


0.38 


— 


- 


-0.10 


— 


— 


7f 


60 3 


14 


3 




0.17 





1.66 


0.45 






-0.08 






7g 


60 3 


9 


3 


12 


0.17 





1.12 


0.12 


0.02 


0.90 


-0.09 


0.39 


0.21 


8a 


60 3 


6 


5 


5 


0.25 





1.22 


0.44 


-0.07 




-0.11 






8b 


60 3 


9 


5 


5 


0.10 





1.20 


0.26 


-0.03 




-0.11 






9a 


60 3 


6 


8 


2.7 


0.27 





1.00 


0.29 


-0.19 




-0.12 






9b 


80 3 


9 


8 


2.7 


0.27 





1.04 


0.36 


-0.08 




-0.09 






9c 


100 3 


14 


8 


8 


0.17 





1.11 


0.44 


-0.01 




-0.06 






9d 


81 3 


14 


8 


40 


0.35 





0.84 


0.07 


0.39 


0.65 


-0.17 


0.22 


0.11 


10a 


60 3 


6 


12 


12 


0.33 





0.88 


0.40 


-0.05 




-0.09 






10b 


60 3 


9 


12 


12 


0.33 





0.90 


0.33 


-0.04 




-0.09 






10c 


80 3 


12 


12 


12 


0.17 





0.90 


0.31 


-0.04 




-0.07 






lOd 


60 3 


9 


12 


60 


0.33 





0.78 


0.06 


0.72 


0.36 


-0.23 


0.03 


0.004 
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Fig. 2. Evolution of spectral energy in different modes 
in a lin-log plot. The inset shows the same, but here the 
ordinate scale is linear. Run 9c. 



In the next subsection we discuss hydrodynamic as- 
pects, i.e. we assume that TZ m is below the critical value 



for dynamo action, which is between 10 and 20. The dis 
cussion of runs with dynamo action follows in Sect. 3.2. 



3.1. Nature of the large scale flow 

In all cases the solutions have a slowly varying large 
scale component, clearly seen in the horizontal averages 
(denoted by overbars) of the velocity. Figure [3] shows 
that u x (z) and u v {z) vary approximately sinusoidally in 
z and are phase shifted relative to each other by about 
90 degrees. This type of mean flow corresponds to a 
Beltrami wave which is helical (tJ • u is here negative; 
uJ = V x u) and approximately force-free (i.e. it has no 
inertia; u 'Vu = 0). Because of weak compressibility, V-w 
is not strictly zero and therefore u z does not need to van- 
ish exactly, but it is nevertheless very close to zero; see 
the dashed lines in Fig. ||. 

The rms velocity of the large scale flow relative to that 
of the full velocity field is given by U^/U. As the value of 
TZ is increased the relative importance of the large scale 
flow diminishes; see Fig. [| This gives us a first indication 
that in the astrophysically interesting limit of TZ — > oo 
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Fig. 3. The three horizontally averaged velocity compo- 
nents for Run 2a (TZ = 1, kf/ki —6) and Run 9c (TZ = 8, 
k{/ki — 14). The relative kinetic helicity of the depicted 
large scale field is = -0.99 for Run 2a and -0.79 for 
Run 9c. 



the AKA-generated large scale flow will be less prominent 
than in the case where TZ is small. 
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Fig. 4. The relative importance of the large scale flow, as 
measured by Uu^jU, versus 1Z. Different values of kf/k\ 
are indicated by different symbols. Note the segregation 
into two different groups. 



Next we analyse the relation between the mean (large 
scale) flow and components of the Reynolds stress, Rij = 



u^u'j , where primes denote deviations from the mean flow, 
i.e. u 1 = u — u. To lowest order the AKA effect couples 
the components of the Reynolds stress linearly to the mean 
flow, 



const 



jkUk + higher order terms, 



(12) 



Table 2. Normalized fit coefficients qo, q\, and q2 for 
models with different values of TZ. 



Run 


fc f /fci 


TZ 


q a /{vk t ) 2 


q± x (vk f ) 


q 2 x (vkf) 


la 


a 



U.oz 


U.ol 


n to 
U. / Z 


U.oo 


2a 


6 


1 


0.25 


0.73 


0.35 


3a 


6 


1.2 


0.35 


0.62 


0.28 


4a 


6 


1.5 


0.52 


0.52 


0.23 


5a 


6 


2 


0.81 


0.41 


0.19 


6a 


6 


2.3 


0.98 


0.34 


0.17 


7a 


6 


3 


1.19 


0.18 


0.07 


7b 


6 


3 


2.23 


0.09 


0.03 


8a 


6 


5 


4.11 


0.043 


0.012 


9a 


6 


8 


6.99 


0.016 


-0.003 


9c 


14 


8 


7.20 


0.009 


0.005 


10a 


6 


12 


10.3 


0.010 


0.002 



which closes the equations for the mean flow. Similarly to 
the alpha effect in dynamo theory, the onjk term produces 
a linear instability for wavenumbers k/kf < ^1Z 2 (Frisch et 
al. 1987). The growth rate is maximum when fc/fcf = jlZ 2 . 
For example, when TZ 2 =2/3 (the value considered by 
Frisch et al. 1987) then the fastest growth is for k/k{ — 
1 /6, which was also the smallest wavenumber ratio in their 
model. Moreover, in this case, this is also the only unstable 
mode; k/kf = 1/3 is already stable and there are no other 
discrete modes in-between. 

For small values of TZ, Frisch et al. (1987) found that 
when the large scale flow reaches saturation, the Reynolds 
stress tensor depends nonlinear ly on u. In particular, its 
horizontal off-diagonal components are given by 



Ra 



in 2 



l + u y /(vkf) + \v? J{vki) 



R 



It j 2 



l-u x /(vki) + \u 2 x /(vk t ) 2 ' 
For small values of |lt| one recovers Eq. ( |l2|) with 

- a yzx = \U 2 /{vk { ) 



(13) 



(14) 



*-xzy 



\1ZU Q . 

In order to check Eqs (|l^ ) and ([lj 
of TZ we generalize this relation to 



(15) 

for different values 



1 



qiUy 



92 u 2 



Ryz - 



Qo 



1 - qiu x 



Q2U 2 X 



(16) 



where go, ?ii and 52 are functions of TZ that can be deter- 
mined by numerically fitting R xz to u y and R yz to —u x . It 
turns out that in these two cases the three fit coefficients 
are the same to a good approximation, in agreement with 
the results by Sulem et al. (1989). We have therefore deter- 
mined the best fit for the combined data set. The results 
are shown in Table || and Fig. S. 

Again, expanding Eq. (p6J) we have 



Rxz ~ go - qoqiuy, 

and therefore 



Ru 



qo + qoqm x , 



a 



xzy 



- a 



yzx 



qoqi- 



(17) 



(18) 
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It turns out that to a good approximation the values of 
QoQi/ivkf) and 2go?2 agree with each other; see Fig. pi 
For TZ < 2, q qi/(vkf) increases with TZ = Uo/{vk{). 
This implies that in dynamical units qoqi is constant: 
Q0Q1/U0 ~ 0.17 for TZ < 2. For larger values of 1Z, 
qQqi/(vk{) decreases like 1Z~ X . In dynamical units, q$qi 
decreases quadratically: qoqi/Uo « 0.61Z~ 2 . 
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Fig. 6. The TZ dependence of qo/(vkf) 2 , qoqi/(vkf), and 
2<3 , o<?2- 




= 3.0 




Fig. 5. i?2; Z vs (plus signs) superimposed with R yz vs 
— u x (diamonds). The solid line gives the fit Eq. |[!) with 
the parameters qo , qi , and 92 from Table |^. 

As we have already seen in Fig. ||, the solutions change 
character near TZ = 3. It turns out that in some cases with 
TZ > 3 the solutions come in the form of travelling waves; 
see Fig. 0. For some values of TZ, however, we also found 
standing waves that alternated in time. A typical value 
for the period is uoT = 600. 

3.2. Magnetic fields 

We now discuss runs where the value of TZ m is large enough 
to allow dynamo action. A good example is Run 7g where 
TZ m = 12. Here the kinematic growth rate of the magnetic 
field is A = 0.02w. The initial field has grown by a bit more 
than five orders of magnitude before it reaches saturation, 
which is when cut w lnl0 5 /0.02 « 600. One clearly sees 
that at this time the large scale flow becomes strongly 




-20 ■ 



3000 



Fig. 7. Time-space diagram of u x (z, t) for Run 7f. fcf/fci 
11. TZ = 3. Positive values of u x are shown in light and 
negative values in dark. Note the migration in the positive 
z-direction. 

suppressed; see Fig. |s|, where we have plotted the evolution 
of kinetic energies contained in various Fourier modes. The 
normalized magnetic field strength, B, saturates at a value 
similar to the normalized velocity, U. Looking at Table 
it is clear that for all runs with a magnetic field the large 
scale flow is strongly suppressed, i.e. U\s,jU never exceeds 
the value 0.1 . . .0.2. 

Prior to saturation of the magnetic field the kinetic 
energy spectrum looks like that in Fig. |], with peaks at 
both the forcing scale as well as the largest scale of the 
system; see Fig. ||. (The multiple peaks to the right of 
the forcing wavenumber come from higher harmonics that 
are here more pronounced than at higher values of TZ; cf. 
Fig. [l].) However, as the magnetic energy approaches sat- 
uration, the kinetic energy becomes suppressed at large 
scales. Figure ^ shows that, in the saturated state, the ki- 
netic energy spectrum has lost its second peak at k = k\. 
Instead, the magnetic energy spectrum has now attained 
a shape similar to that of the kinetic energy before satu- 
ration. In that sense, the inverse cascade-type behaviour 
in velocity is now replaced by an inverse cascade-type be- 
haviour of the magnetic field. 
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Fig. 8. Evolution of spectral energy in different modes 
in a lin-log plot. The inset shows the same, but here the 
ordinate scale is linear. Run 7g. 
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Fig. 9. Power spectra of kinetic (upper panel) and mag- 
netic (lower panel) energies for Run 7g at different times 
before saturation of the large scale fields (dotted lines) and 
after saturation (solid lines). Note that after saturation 
the inverse cascade-type behaviour in velocity is replaced 
by an inverse cascade-type behaviour of the magnetic field. 



In Fig. [Oj we have plotted the evolution of U, Uls, 
and B (upper panel) together with the relative kinetic and 
magnetic helicities, 7Yk and Hmj as well as the nondimen- 
sional length scale related to the magnetic helicity, 7Ym 
(lower panel). The kinetic helicity is always negative and 
temporarily enhanced around the time when the magnetic 
field reaches saturation. We recall that at each instant the 




1500 



Fig. 10. Evolution of U, Uls, and B (upper panel) and 
T~Ck, T~Cm, and Hm (lower panel). Run 7g. 

forcing has zero helicity, and it is only due to the temporal 
shift of the forcing pattern that helicity is introduced into 
the flow. 

The magnetic field has positive helicity (Fig. [h]). The 
opposite signs of kinetic and magnetic helicities are in 
agreement with what is found for helically forced turbu- 
lence; see Brandenburg (2001), where the forcing was how- 
ever chosen to have positive helicity. This led to positive 
helicity of the flow, so all signs are reversed compared to 
the present case. Indeed, closer inspection of the power 
spectrum of the magnetic helicity shows that at small 
scales the signs of kinetic and magnetic helicities agree. 
However, because of (approximate) helicity conservation 
the integrated helicity spectrum has to vanish, which is 
why at large scales the sign of magnetic helicity is re- 
versed. 

4. Conclusions 

The present calculations have verified the possibility of the 
formation of large scale flows from small scale flows lacking 
parity invariance. These large scale flows tend to be helical 
and of Beltrami type, as expected from the nature of the 
anisotropic kinetic alpha (or AKA) effect. The power of 
the large scale flows tends to be more strongly distributed 
over several Fourier modes. This is in stark contrast to 
the magnetic inverse cascade or the dynamo alpha effect 
which tends to select just one large scale Fourier mode - at 
least in the case of a fully periodic box; see Brandenburg 
(2001). 

The resulting flows support the possibility of dynamo 
action once the magnetic Reynolds number exceeds a value 
typical of non- helical dynamo action. This critical value of 
the magnetic Reynolds number, based on the wavenumber 
of the forcing, is around 10. For comparison, in the case 
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of fully helical forcing, the critical dynamo numbers are 
about ten times smaller; cf. Brandenburg (2001), where 
the magnetic Reynolds numbers based on the forcing scale 
need to be divided by 2tt to be compatible with the present 
work. However, once the magnetic field reaches saturation, 
it begins to suppress the large scale flow. 

In the present work we have adopted the forcing func- 
tion of Frisch et al. (1987). This has the advantage that we 
were able to make contact with earlier studies where the 
AKA effect was also considered numerically. This flow was 
constructed on purely mathematical grounds in order to 
demonstrate the very existence of the AKA effect. Other 
more realistic types of forcing have been proposed by 
Levina et al. (2000). It is important to extend the present 
studies to these types of forcing in order to sec whether the 
AKA effect works for broader classes of forcing. However, 
our present results suggest that the anisotropic kinetic al- 
pha effect should not play a role in astrophysics where the 
value of 1Z is in general large and magnetic fields are usu- 
ally present. It should be emphasized that in some of the 
best cases where a large scale velocity pattern emerged 
(e.g. Run 7f), the large scale flow pattern is never really 
as pronounced as the large scale magnetic fields that are 
produced by helical turbulence (cf. Brandenburg 2001). 

Finally, we should mention that our results do by no 
means address the question of whether or not large scale 
vortices can form astrophysical bodies such as giant plan- 
ets or accretion discs. Such vortices are long-lived, quasi- 
stable formations, possibly belonging to the class of so- 
lutions studied by Goodman et al. (1987). It it possible 
that they are formed simply as a matter of suitable ini- 
tial conditions, but they were also found in simulations of 
Hawley (1987). This type of solution would be essentially 
nonlinear, in contrast to solutions of the mean-field equa- 
tions with AKA effect which are possible already in linear 
approximation. 
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